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Polishing the Lens: I Pionic Final State Interactions and HBT Correlations- 
Distorted Wave Emission Function (DWEF) Formalism and Examples 

Gerald A. Miller and John G. Cramer 
Department of Physics, University of Washington 
Seattle, WA 98195-1560 

The emission of pions produced within a dense, strongly-interacting system of matter in the 
presence of strong radial flow and absorption is described using a relativistic optical model formalism, 
replacing the attenuated or unattenuated plane waves of earlier emission function approaches with 
"distorted wave" solutions to a relativistic wave equation including a complex optical potential. The 
resulting distorted-wave emission function model (DWEF) is used in numerical calculations to flt 
HBT correlations and the resonance-corrected pion spectrum from central-collision STAR Au-|-Au 
pion data at y/s = 200 GeV. The parameters of the emission function are constrained by adopting a 
pion formation temperature taken from lattice gauge calculations and a chemical potential equal to 
the pion mass suggested by chiral symmetry restoration. Excellent agreement with the STAR data 
are obtained. The applications are extended by applying linear participant scaling to the space-time 
parameters of the model while keeping all other parameters fixed. This allows us to predict HBT 
radii over a range of centralities for both Au-I-Au and Cu-|-Cu collisions. Good agreement is found 
with STAR HBT data for all but the most peripheral collisions. The squares of pionic distorted wave 
functions, obtained as exact numerical solutions to the wave equation, are displayed and significant 
differences with the results of using the familiar eikonal approximation are found. Using the eikonal 
approximation leads to a qualitative accounting for the effects of the imaginary part of the optical 
potential, but fails entirely to include the effects of the real part of the optical potential. A simple 
example is used to illustrate that an attractive optical potential can have large effects on extracting 
radii and can also lead to oscillations in radii measured at low momenta. 



I. INTRODUCTION 

Measurements of the two-particle momentum correlations between pairs of identical particles have been 
used to study the space-time structure of the "fireball" produced in the collision between two heavy ions moving 
relativistically. The quantum statistical effects of symmetrization cause an enhancement of the two-boson 
coincidence rate at small momentum differences that can be related to the space-time extent of the particle 
source. This method, called HBT interferometry, has been applied extensively in recent experiments at the 
Relativistic Heavy Ion Colhdcr (RHIC) by the STAR and PHENIX collaborations. See the reviews[l|l|l|l|. 

The invariant ratio of the cross section for the production of two pions of momenta pj^, p2 to the product of 
single particle production cross sections is analyzed as the correlation function C(pi, P2). We define q=Pi-p2 
and K=(pi-|-p2)/2, with Kt as the component perpendicular to the beam direction. (We focus on mid- 
rapidity data, where K = Ky-) The correlation function can be parameterized for small q as C(q, K) — 1 w 
Aexp {—RoQo ~ ^sls ~ ^lIl) ~ ^ ^olo ~ ^sls ~ ^lIl) (liRi ^ !)> where O, S, L represent directions 
parallel to Ky, perpendicular to K.t and the beam direction, and parallel to the beam direction^. Early and 
recent Q hydrodynamic calculations predicted that a fireball evolving through a quark-gluon-hadronic phase 
transitions would emit pions over a long time period, causing a large ratio Rq/Rs- The puzzling experimental 
resuh that Rq /Rs « 1 is part of what has been called "the RHIC HBT puzzle" |i| . 

Another part of the puzzle is that the measured radii depend strongly on the average momentum K, 
typically decreasing in size by about 50% over the measured range. This dependence of a geometrical parameter 
on the probe momentum shows immediately that the radii are not simply a property of a static source. The 
influence of the interactions between the pionic probe and the medium, as well as the effects of transverse flow, 
must be taken into account when extracting the radii. The medium at RHIC seems to be a very high density, 
strongly interacting plasma[3 , so that any pions made in its interior would be expected to interact strongly. 

We studied the effects of including the pionic interactions in previous workjlOf, finding that the only way 
to simultaneously describe the measured HBT radii and pionic spectra is to include the effects of pion-medium 
final state interactions by solving the relevant relativistic wave equation. These interactions are so strongly 
attractive that the pions can be taken as propagating through a system with a restored chiral symmetry. 
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The principal aim of the present work is to present a detailed treatment of the formalism that will allow 
wide application of our technique, and we present some new applications here. We also provide specific simple 
examples to demonstrate that the effects of pionic interactions cause the measured sizes of the medium to 
be different than the true sizes. Furthermore, we shall explicitly demonstrate that classical treatments of the 
pion- medium interactions, based on using the cikonal approximation for solutions of the wave equation arc not 
valid. 

An outline of the remainder of this paper follows. Previous standard formalisms that use plane wave pions 
are briefly reviewed in Sect.^] The technical method of incorporating the influence of final state interactions 
between the pion and the medium is described in Sect. IIIII Pionic emission in the absence of final state 
interaction is described with an emission function So that is of a form motivated by hydrodynamics. The 
connection between the symmetries of this function (for the case of head-on collisions) and the form of the 
pion optical potential is described in Sect. IIVI The role of the complex optical potential and the use of chiral 
symmetry to constrain its form at low energies is explained in Sect. ^ The the specific numerical algorithm 
necessary to incorporate the optical potential is presented in Sect. IVII Once our approach is defined, there 
is only one approximation we need to make. The validity of this requires only that the source size be much 
larger than the inverse of the temperature. This large source approximation, LSA, is explained in Sect. IVIII 
The resulting distorted wave emission function, DWEF, is evaluated using two different equivalent methods in 
Sec. IVIIII In Sect. |IX| wc apply these methods to STAR central Au+Au data at a/s = 200 GeV. Instead of 
treating the temperature of the system as a free parameter, as was done in |lOl| . we fix its value at the transition 
temperature = 193 MeV, obtained in the most recent lattice QCD calculation [Tll|. 

For conventional calculations of the spectra, chemical equilibrium analyzes yield lower temperatures Tch = 
174 MeV ^ large difference between Tc and Tch implies that the hadrons interact after the dcconfinement 
transition occurs. This notion is entirely consistent with our treatment of pionic distortions which has as its 
fundamental assumption that pions interact in a hot dense medium before escaping to freedom. We also fix 
the pion chemical potential at the pion mass (/z^ = 139.57 MeV). The resulting procedure reduces the number 
of free parameters by two to a total of nine. Then we are able to reproduce a total of 32 data points with 
high accuracy. Section Hxl also extends earlier numerical results^3 by computing the dependence of the HBT 
radii on the centrality in Au+Au collisions and by making predictions for Cu+Cu collisions vs. centrality. The 
eikonal approximation to solving the wave equation is discussed in Sect. Q in which the important influence of 
the opacity and the vanishing of the effects of the real potential arc also described. The importance of the real 
part of the optical potential in obtaining oscillating radii at low energies is illustrated through two examples in 
Sect. IXl| A brief summary is presented in Sect. IXIll Finally, a short appendix verifies an approximation to an 
integral. 



The aim is to include the effects of final state interactions of outgoing pions. We'll begin with a brief 
review of the formalism previously used to describe HBT correlations for situations in which the pions do not 
interact with the medium. 

The relevant observables are the covariant single- and two-particle emission functions defined as appro- 
priately normalized ratios of cross sections [T^ 



II. PREVIOUS FORMALISM - PLANE WAVE PIONS 



V2{P1,P2) ^ E1E2 




(1) 



(2) 




(3) 



(4) 



C(Pl,P2) 



(iV)2 7'2(Pl,P2) 7'2(Pl,P2) 



(5) 



{N{N - 1)) 7'i(pi)7'i(p2) Pi(pi)Pi(P2) ■ 



The last step of Eq. ||SJ) is obtained because {N'^) S> {N) for the very high energy collisions of interest here. 

These observables can be expressed in terms of an emission function S{x, K) that is the Wigner transform 
of the density matrix associated with the currents that emit the pions. The literature [^Il4j presents the single- 
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particle emission function as 



So{x,p)= I ^^exp{~ip-y){r{x + y/2)J{x-y/2)) 



2{2ny 

where p is the four-momentum, of an on-shell emitted pion and J{x) is the current that acts as a source of pion 
fields. We shall assume that J{x) is not altered by the presence of the produced pions. The brackets indicate 
that one takes an ensemble average over chaotic sources. Thus we assume that the emission process is initially 
uncorrclated: the pions are emitted from chaotic sources that have random phases [l^ IT^ ITsI l . This assumption 
seems consistent with analysis of data produced at SPS and RHIC0 The second term of Eq. © is written 
to show that the emission function is determined by the product of a plane wave factor and the current that 
emits the pions. This is the standard result of scattering theory, if final state interactions are ignored. Here 
and throughout the paper we use natural units in which h and c are unity. Our use of the subscript denotes 
that the effects of final state interactions arc ignored. 

Next we discuss the relationship between the function So{x,p) and the single-pion production cross 
section. If the source current couples linearly with the pion field (with a chosen interaction Lagrangian, Lj), 
and no final state interactions occur, the matrix element of the interaction —Lj between the vacuum and 
single-pion final state is 

J{p) = J d^xe'P '^Jix) = J d'^x(j)p{x)J{x), (7) 

where (j)p = e^P'^ is the solution of 

^-V^+m^)0p = O. (8) 
In this plane wave (pw) approximation, the covariant single-particle emission function can be expressed as 



vr{p)^{\j{p)n = J d\so{x,p), (9) 

The literature 01 presents the two-particle emission function as 

Soix, K)= J exp i-iK ■ y){J*{x + y/2)J{x - y/2)). (10) 

The quantity that enters in the HBT correlation is 5*0(2;, q) with 

So{x,K,q) ^ So{x,K)e-^^-'^ - J -^(e-f-(-+«/2) J*(.x + y/2)e^P-(--^/2) J(x - y/2)> ^ 5o(x,pi,P2), (H) 
where pi,P2 are the four-momenta of the two on-shcU emitted pions, and 

q = Pi-P2, K = ^{pi+P2)- (12) 

The brackets again indicate that one takes an ensemble average over chaotic sources. The second form, Eq. (Illll 
shows again that products of plane-wave functions and currents determine the observables. The notation 
Sq{x,Pi,P2) is introduced to emphasize the two-particle nature of the emission function. 
In plane wave approximation the correlation function is given by the expression 



vriPi)vriP2)c,ipuP2) = ( 



/ |g3;^(e- -e--+e--e--) J(.,)J(x.)' 



),(13) 



in which the effects of the Bose statistics is explicit. One proceeds by taking the ab- 
solute square to find that the expression contains four terms. The direct terms involve 

(^e*Pi-^ie-'Pi-^'i(J(xi)J*(a;'i))e*^'2-^2g-v2-4(j(a;2)J*(a;2)) plus 1 ^ 2^ and the exchange terms involve 
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(^e*Pi-^ie-*P2-^'i(J(a;i)J*(x'i)))e*J'^-^^e-*Pi-^2(J(a;2)J*(x2))) plus 1 ^ 2^ . The placement of the brackets 

arises from the use of chaotic sources: the effects of pions produced by two different chaotic sources average to 
o[il[i|. Using the absolute square and the previous definitions in Eq. (|13|) yields the result 

Cq{Pi,P2) = I + ^ ,2 , ~ (14) 



III. INTRODUCING FINAL STATE INTERACTIONS 



Our emission function, 5*0, is not meant to be the same as that obtained in conventional hydrodynaniic 
calculations of hadronic spectral functions. In such calculationsQ the quantity 5*0(2;, A') is assumed to be a 
local equilibrium Bose-Einstein distribution localized on a 3-dimensional frcczc-out hypcrsurface that separates 
the thermalized interior of the hot dense medium from the free-streaming particles on its exterior |20l|. Here 
we assume that the emitted pions (and other hadrons) undergo significant interactions while escaping the hot 
dense medium. The emitted pions interact both with the hot dense matter and other hadrons as they escape the 
system. This is emission from a coexistence phase. The purpose of the present section is to provide a formalism 
that allows such effects to take place, while also including the usual effects of emission from a freeze out surface. 



A. Distorted Waves 



We wish to include the effect that an escapi ng p ion has to "fight" its way through the medium. Our 
formalism relies heavily on earlier derivations in Ref. |14| . The effects of the fighting distort the wave away from 
its plane wave form dictated by Eq. ijHJ). These interactions of the pions with the medium, represented by the 
optical potential [/, lead to a modified equation (approximated as a one-body equation) 

(^-V^ + LZ + m^^ V = 0. (15) 

This equation provides an approximate treatment of the complicated final state interactions. In principle, 
the optical potential should be related to the underlying pion source currents. We adopt a phenomenological 
approach. The resulting S-matrix depends on the exact solution, ipp \x)^ to this equation with the out boundary 
condition of approaching at t ^ 00, the free wave ^p(a;)|23|. 

Ref. [l4| uses standard S-matrix reduction techniques to show that matrix element of the pion field 

between the vacuum and a one-pion final state p is simply the complex conjugate of ipp \x)^ if the optical 
potential is introduced as in Eq. H15|l . This means that in calculating the matrix elements of — i/ between the 
vacuum and multi-pion final states one simply replaces the plane wave solutions of Eq. (jSJ by the out-going 
wave solutions of Eq. H15() . Thus the simple replacement: 

e'P-^^(,-)*(.^) (16) 

is sufficient to include final state interaction effects within the optical model approximation. Therefore the 
single-particle emission function is given by 

S{x,p,) = J i^{r{^ + ymJ{^ y/2))^^;^(x + y/2)47)*(2; - y/2). (17) 

This expression reduces to the plane wave form Sq{x^p\) Eq. H10() if the distorted wave i/jpj^ are replaced by 
plane waves. The two-particle emission function that includes the effects of final state interactions is obtained 
by replacing Eq. H1U|I by the result 

S(x,p2.V^)^S{x,K,q)^ J -^(J*(x + y/2)J(x--2//2))V'(;)(^ + y/2)4;'*(x--2//2), (18) 



which applies for calculating the two-particle emission function. In the plane wave limit, S{x, K, q) 
So{x,K)e-"i-''. 
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Physical observables for the emission of two pions of momenta pi , p2 are determined by the correlation 
fmiction C{q, K) [K = ^{pi + P2), q ^ pi — P2), which is given by 

^'^' ^ ^ S d^xS{x,pi) J d^xS{x,p2y ^ ' 

One obtains the usual expression (|14|l if the wave functions "i/jp ^ are replaced by plane wave functions ('0p ^ {x) — > 

The expression (|18|) contains the ensemble average of the currents. This may be expressed in terms of 
Sq{x^ K) by taking the Fourier transform of Eq. ©. We then obtain the convolution formula: 

S{x,K,q) = j d*K'S„{x,K') I -0^ e^'^'-y^l-^x + y/2)^(^->{x ^ y/2), (20) 

where the subscripts indicate the momenta pi,P2 of the detected pions and K = ^{pi +^2)- The quantity 
S{x, K, q) is used to compute experimental observables in the same way that was previously used. This is 
the distorted wave emission function DWEF formalism. The expression H2U|) is the coordinate space version of 
Eq. (5.25) of Ref. [H. 

We emphasize that Sq (which reflects the true properties of the source) is no longer directly related 
to observables-the appearance of distorted waves obscures the relationship between the data and the true 
properties of the source. 



B. Extracting HBT radii 

The correlation function of Eq. are related to HBT radii in two different ways. In the first method, 
one treats the momentum differences qo,s,L as small quantities and then expands keeping terms to second order 
so that: 

c{q,K)-i^i-q'oRl-QlRl-qlRl, (21) 

where qo is the transverse component that is parallel to the direction of K, qg is the transverse component that 
is perpendicular to the direction of K, and qL is the longitudinal component. Here and below, because we are 
focusing on central collisions we ignore the qoQLRoL cross term. Another parameterization is 

C{q, K) - 1 « c^vi-qlRl - 4RI - qlRl), (22) 
In practice, describing data and extracting radii require using: 

c{q,K) - 1 « XcM-q^Rl - ilRl - qlRl), (23) 

C(g, A') - 1 « A(l - q'oRh - 4RI - llRl)- (24) 

Here the reduction factor A (typically about 1/2) is the fraction of pairs that originate in the space-time region 
relevant for correlations, see the review 0. Current understanding is that the HBT data are consistent 
with incoherent emission, and accounting for the many pions that are produced by the decays of resonances 
far outside the collision region can reproduce the factor A. These "halo" pions do not have a BE-enhanced 
correlation in the q region measured with the pions emitted from the core of hot dense matter, but they cannot 
be experimentally separated from the latter. 

Our calculations employ the core- halo model |l(il | in which pions are assumed to arise either from the hot, 
dense core or from the halo. This model is a simplified version of more detailed treatments of resonance decays 
that are discussed in the review . 

The infiuence of resonances affects the extraction of radii and the measurements of the pion spectrum 
in different ways. We account for this in our phenomenological analysis, see the erratum of Ref. [lol | and Sect. 
IX-A. Here we note only that effects of pions produced by short-lived resonances (such as A, p, • • •) that are not 
explicitly included in the pionic wave equation are included in the function 5*0 . The pions resulting from long 
lived resonances are not included in Sq. Pions resulting from the decay of slowly moving U, mesons that occur 
inside the dense matter system are included in 5*0, but the pions from rapidly moving lo mesons that decay 
outside of the system arc not included in Sq. 
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The approximate forms (|23I24|I provide two ways to extract radii. The former suggests that 



1 

C{quK) 



= —^^TuZ-i7\ 7' i = 0,S,L (25) 



while the latter implies 



i^^-4 ^^^-f "\ -0,^,L (26) 

Eq. (|25|l can be used for values of qi such that (C — 1)/A can be approximated by a Gaussian function. The use 
of Eq. H24|l requires that qiRi ^ 1. We show below in Sect. IVIlTl that. while the correlation functions are not 
Gaussians (so that the squared radii are not moments of the correlation function), the Gaussian parameterization 
is quite accurate at the modest (but not very small) values of qi ~ 30 MeV/c that dominate the experimental 
extraction of radii. 



IV. SYMMETRIES OF So{x,K) AND THE FORM OF THE PION DISTORTED WAVES 



We use the hydrodynamic parameterization of the source of Ref. 0, ITsl ITflj l . Generally based on the 
Bjorkcn tube model, it is given by 



S'o(x, K) ct^x 



Mj_ cosh(77 - Y) 
(2^ 



Xrdr 



VMM' 

rdr drj b db d(j). 



exp 
exp 



- 1 



p{b) exp 



2(A77)2 



(r-roY 



2(At)2 



di] b db ( 



(27) 
(28) 



b are equivalent to x± , 77 = i In , r 



Here, /.t^ is the pion chemical potential, the variables (6, 4> 
\/W~~z^, Mt = V K\ + m\ , F = i In g^^^' . The factor in the brackets involving r has the same normalization 
as the delta function: (5(t — tq). We use a Bose-Einstein distribution instead of the Boltzmann distribution 
of Ref. ^3 and also allow the transverse density p(&) to have a general form instead of a Gaussian. Note 
that the model contains the parameters: i?, 770, A77, Ar, tq, 77/. As stated originally, temperature gradients were 
included, but we treat the temperature as a constant for our numerical calculations. We shall see below that 
our formalism can easily be generalized to allow the temperature to vary as a function of 6. 

The longitudinal 77, r and transverse (b, K) variables can be separated in the following form: 



5*0 (x, if) 
B,(b,K) 



5o(77,T,r) B^(b,K) 



exp[(/v • u - ii^)/T] - 1 



P{b) 



Q I v\ - cosh(77 - Y) 



(27r)3 



exp 



exp 



2(At)2 



(29) 
(30) 

(31) 



2(A77)2J ^2^(Ar)2 

Here p{b) is a function, normalized as p(0) ~ 1 that represents the transverse density. To be specific we use 

p{b) = [l/(cxp((6 - Rws)/aws) + !)]'• (32) 

This distribution has a correct exponential fall-off at large values of fe, and different choices of the parameters 
Rws, CLws allow a variety of shapes to be assumed. 

We concentrate on the kinematics of the STAR experiment which detects pions within half a unit of 
rapidity of moving perpendicular to the beam, so we take y = 0. This means that the average momentum is 
transverse: K = Kj_ = Kt. With colliding beams of equal mass and energy there is a fore-aft symmetry along 
the longitudinal axis, so that we use 7/0 ~ 0. Therefore it is useful to define 



So{il,T) = So{ri,T,Y = 0). 



(33) 
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The velocity field u(x) describing the dynamics of the expanding source is parameterized by |22 | 

u^{x) = (cosh ?7 cosh ryt(&), cosi/i sinh?7t(6), sin(/> sinh?7t(6), sinhrycoshyyt (&)) , (34) 

with (j) the angle between and u (or when appearing in the single-particle emission function, it is the angle 
between and u). Eq. H34|l implements a boost-invariant longitudinal flow profile vl = z/t, with a linear 
radial profile of strength r]f for the transverse flow rapidity: 

m{h)=Vf^- (35) 
The exponent that enters in the Bose-Einstein distribution is given by the four- vector dot product: 

K ■ u{x) = AI±^ cosh(r/ — Y) cosh rjt (b) — K± sinh i]t (b) cos (j) . (36) 

The presence of a non- vanishing value of 77/ causes the emission function to depend on both K and M±. The 
Bose-Einstein distribution of Eq. (|3()|l is evaluated as a sum of Boltzmann distributions: 

, ' , = E exp( -^"; + ^- ), T„ . I. (37) 



exp 



T 



The parameterization H27() is motivated by hydrodynamical models with approximately boost-invariant 
longitudinal dynamics. It uses thermodynamic and hydrodynamic parameters and appropriate coordinates. 
The "emission function" given above was originally intended to parameterize the distribution of points of last 
interaction in the source. In conventional treatments, the CooDer-Frve|20j matching procedure is used to obtain 
distributions of detected particles. Our approach is more general. We assume that pions can be formed at 
any space-time point (including but not limited to the freeze-out surface) during the collision, and propagate 
through the dense medium while interacting before being detected. Thus we take S'o(a;, A') to be the emission 
function in the absence of final state interactions. 

We note that the emission function of Eq. (|27|l has been previously used in the "blast- wave model" , 
and we would like to comment on the differences between our formalism and that model. These include: 

• The blast-wave model does not attempt to reproduce the normalization of the spectrum, and sets the 
chemical potential [i^; to 0. We find that taking [ij, to be the pion mass allows us to reproduce both the 
normalization and the shape of the pion momentum spectrum. 

• The blast-wave model uses the smoothness approximation and computes HBT radii as moments of the 
emission function. 

• The blast-wave model uses plane waves and therefore omits the effects of the optical potential. 

• The blast-wave model uses a parameterized emission function to describe the six-dimensional distribution 
of pions on some freeze-out hypersurface, while the emission function of the DWEF model describes the 
initial emission of pions within the hyper-volume of the hot, dense medium produced by the collision. 
Therefore, it is not appropriate to compare emission-function parameters like T and R that are derived 
from blast-wave fits with those of the present model. 



V. FINAL STATE INTERACTIONS AND THE OPTICAL POTENTIAL 



The salient feature of the 200 GeV data is the high density of the produced matter, so we treat the 
effects of pion final state interactions with a dense medium. We adopt a single-channel approach that uses the 
interaction-distorted incoming wave '^p^^* (xi). 

We assume that the matter produced in the central region of the collision is cylindrically symmetric with 
a very long axis, so that an expression of the form 129|) is valid. In that case, the optical potential U representing 
the interaction between a pion and the medium is a complex, azimuthally-symmetric function depending on 
pion momentum and local density. Within our formalism the influence of some time-dependent effects in U 
introduced by the time-dependent source So is incorporated in the energy dependence of the optical potential, 
and the pion- medium interaction time is restricted by Sq. 
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The optical potential accounts for the interaction between each pion and the surrounding medium, but 
does not include the interaction between the two pions. In particular, the Coulomb interaction is known to be 
important. The experimental analysis removes this effect before extracting radii from the data. The Coulomb 
interaction between pions is of long range and its important effects occur when the pions are outside the medium. 
Indeed, specific analyzes show that Coulomb effect occur only at very low relative momenta In contrast, 
the optical potential acts only when pions are inside. We therefore expect that presence of the optical potential 
would not influence the removal of the Coulomb interaction. A quantitative treatment of the effect of the optical 
potential on the removal of Coulomb effects would involve solving a three-body problem. This is not warranted 
at the present stage of development, but might eventually become worthwhile. 

The optical potential accounts for situations in which the pion changes energy or disappears entirely due 
to its interactions with the dense medium. We do not assume to know the content of the dense medium, and 
therefore will use a phcnomcnological optical potential. But it is worthwhile to consider a simple example to 
get an idea about how large the optical potential can be. Suppose, e.g., that the medium is a gas of pions. 
Then tttt scattering would be the origin of U. In the impulse approximation, the central optical potential would 
be Uo = — 47r/po, where / is the complex forward scattering amplitude and po the central density. For low 
energy pion-pion interactions, 47rlm[/(p)] = pa, with cr « 1 mb. At a momentum p = 1 fm^^ = 197.3 MeV/c, 
using a pion density about ten times the baryon density of ordinary nuclear matter. Im[C/(0)] ~ —0.15 fm~^, 
representing significant opacity. 

The optical potential must be an analytic function of energy, and therefore the existence of an imaginary 
part mandates the existence of a real part. Thus, any analysis needs to treat U as a complex function. Under 
certain circumstances the real part can be very large. For example, if two interacting pions each have less energy 
than half of the rho meson mass, the final state interactions caused by virtual transitions to a rho meson would 
be strongly attractive. Additionally, the influence of chiral symmetry restoration can lead to a strong real part. 
This is discussed next. 



A. Chiral Symmetry Restoration 

Suppose the dense medium is one in which chiral symmetry is restored. This means that the value of the 
quark condensate vanishes, an effect that could be caused by an increase in temperature or density. The pion 
mass is proportional to the quark condensate via the GMOR relation [23| 

mlf^ = -^^^^^±^{0\uu + dd\0), (38) 

where /jr is the weak pion decay constant w 93 MeV. It is believed that m^, /^r and the condensate {0\uu+dd\Q) all 
depend on temperature and density. If one takes the Brown-Rho|22 scaling relation for fj^ and the perturbatively 
calculated temperature dependence of the condensate, the pion mass is proportional to the cube root of the 
condensate, and therefore vanishes for sufficiently large temperatures. See the reviews l^^- Suppose the optical 
potential arises only from the temperature dependence of the pion mass. Then the Klein-Gordon equation 
would take the form; 

(-V2 + m2(r))V;-(p2 + m2)V>, (39) 

for regions inside the medium. The effects of the medium are incorporated through the difference between 
m — 7r^(T) and mj. If one rc-writes Eq. (p??^ as a Klein-Gordon equation (— + U)iIj = p^tp, then the optical 
potential takes the form: 

Uib) = {ml{T)^ml)p{b), (40) 

in which the finite extent of the medium is accounted for by the factor p{b). If mT^[T) approaches zero, the 
optical potential is attractive with magnitude mt. 

A more recent study by Son & Stephcnov[33 provides a more detailed treatment of the effects of chiral 
restoration in which the general p-wavc nature of the low-energy interaction between pions and any target is 
included. We are guided by this work. Son & Stcphcnov use the dispersion relation for low momentum pions 
in infinite nuclear matter|30ll3lj: 

LO^ ^ul[f +m^{Tf), (41) 

where fP' is the infinite-sized matter version of — V^. The quantity u^r is termed the pion velocity, even though it 
is only that when m.^{T) vanishes. The term 'mT^{T) is denoted the pion screening mass. This quantity appears 
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in the expression for the static Euchdcan pion correlator y^- The energy of a pion at p = is termed the pion 
pole mass. The free pion mass is 7717^. In Ref. Eq- l|41|l applies only for T < Tc- For larger temperatures, 
chiral symmetry is restored and pions arc massless. 

Defining tEE {Tc- T)/Tc, SS find ml{T) - t^-'', - t^, with P < 1^, eg ly = .73,(3= .38^. These 
equations arc valid for temperature close to (but not too close to) the critical point. Another view about the 
dispersion relation can be found in Ref. |33|| . This general discussion about the influence of chiral restoration 
provides some guidance, but does not tell us exactly to use. 

We wish to obtain an equivalent optical potential and see if it is attractive or repulsive. Use the Klein 
Gordon equation in the form 

+ U + ml =p^ +ml. (42) 

In regions outside the dense mediimr where [/ = 0, the operator is simply the square of the momentum, p^ . 
Subtract from to obtain 

U = ulmliT) -ml + {ul - l)f , (43) 

an expression that is the sum of two negative definite terms. This form can be simplified by using the wave 
equation H42|l to remove term fp' . Then one finds a momentum-dependent optical potential: 

^ ^ ulmljT) - ml + {ul - l)p^ ^ ^^^^ 
ul 

Note that if becomes really small the optical potential becomes very strongly attractive. 

For matter of finite size, the term p^ can also be interpreted as —V • cV which is the Kisslinger term p^ . 
For infinite nuclear matter only forward scattering occurs and the two terms are identical, but differences may 
arise for scattering from media of finite size. One can not tell the difference between the two terms at the start, 
so that we find a general form 

U{b) = -{wo + W2{1- e)p^)p{b) - ew2^p{h) ■ V, (45) 

with both the real and imaginary parts of u'o,W2,e positive for attractive interactions, and p{h) taken from 
Eq. (|32|l . This simple form is strictly valid only for low-energy pions. The limit of an infinite tube is used to 
write U as independent of z. We find for the 200 GeV data, that for temperatures below our assumed value of 
193 McV the fitting prefers a very small value of the Kisslinger term, so we simply set e = 0. We also take wq 
real (so that there is no opacity at p=0). We shall see below, that if the temperature is set to a value greater 
than 193 MeV, the fit is improved by including the gradient terms at about the 20% level (e « 0.2). The simple 
form Eq. H45|l is sufficient to account for the data we study, Kt < 600 MeV/c, but we remind the reader that 
the interaction strength does not grow as p^ for p much greater than about 400 MeV/c. 

We shall discuss the precise parameters of the optical potential in subsequent sections. For now we may 
simply proceed using the assumption that J7 is a complex, azimuthally-symmetric function depending on pion 
momentum and local density of the form of Eq. H45|) . 



VI. FINDING V'*'(x) 

The evaluation of the emission function (|20|) requires performing an eight dimensional integral using a 
distorted wave. We shall use symmetries to reduce the number of numerical evaluations. Doing this depends 
on obtaining a compact expression for the distorted wave. In the present section, we show how the distorted 
waves are evaluated. 

The first step is to realize that the function ipp •* {x) represents an energv-eigenfunction|l4| provided the 
optical potential does not change with time. We shall show below that the value of 5t is not large. Thus the 
the time-independent optical potential that we use can be thought of as a time-averaged optical potential. So 
we have 

V;(,-)(x) = e"''^" ^V(,-)(x). (46) 

To proceed we need to examine the properties of the wave function ?/;p~-'(x). It is conventional to compute 
i/'p^'' (x) , and we will follow this convention. Then we use time reversal invariancc in the form 

vI/(-)(x) = *L+'*(x) (47) 
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to obtain the desired wave function. 

The next step is to reahze that for central colhsions, Y ~ 0, the emission function has a cyUndrical 

symmetry. This means that the expected optical potential is azimuthally symmetric. If we take the matter to 
have the form of a very long tube, the optical potential will be independent of z. Then one obtains a solution 
that takes a product form 

vl/(;)(x) = e™-^/V^;,Ux^=b), (48) 

Pi,2 = K±q/2±zqL/2, (49) 

where the vector q is defined as a transverse vector q • z = 0. 

We may obtain the wave function ipp'^ (t>) by solving the wave equation 

(-Vi + C/(5))4-)*(b) =p2^.(-)*(b). (50) 

If [/ = 0, Vp^^(b) ^ e'P ^- Many previous treatments of opacity can be understood as using the eikonal 
approximation to obtain solutions to Eq. H5U|) . 

We take the optical potential to have the azimuthally-symmetric form of Eq. H45|l so that the solution 
for Vp'^^(b) can be expanded in partial wave form in plane polar coordinates (b -^/x^, (f>, coscj) = p • b): 



4+'(b)= E .^n(p,^)^"e™^ (51) 

m— — oo,oo 

4+' (b) = fo{p, b) + 2 frniP, b)i"' COS m0, (52) 

m — l.oo 

with (|52|l taking into account the invariance of the differential equation for fm under the interchange b —> —6. 
Note that we may use Eq. H47|) to find 

V'^-'*(xi) = /o(p,6) + 2 J2 /n^(p,6)(-^)" cosm<^. (53) 

m— l,oo 

In practice a finite number of terms is needed, with m < rrimax ~ 2pR/h. 
Use Eq. (O in 1^01) to find 

1 d \ 62 

^ + l^, + ip'- -^))UpM - Ui-,)fAP,b) = 0. (54) 

Note that for large enough b, U vanishes and the fm are linear combinations of Bessel Jm and Neumann Nm 
functions or Hankel functions Hm '^\ 

This differential equation can be solved numerically using the Runge-Kutta technique. One determines 
the wave function by matching the numerical solutions to the analytic solution 

/,„ (p, b) = Am ( J™ [pb] + TmH^i'^ (p6)) . (55) 

One matches the numerical function and its derivative (at large enough b so that U = 0) so as to determine 
the constants Am, T„i- The normalization of V'p^'' is such that A„i = 1-asymptotically the wave is a sum of an 
ordinary plane wave and an outgoing wave, and the functions /m(p, b) can be regarded as phase shifted Bessel 
functions. Using the partial- wave form of two-dimensional wave function (|52|l will simplify the evaluation of 
S{x,K) of Eq. (EOJ. 

VII. THE DISTORTED WAVE EMISSION FUNCTION (DWEF) AND THE LARGE SOURCE 

APPROXIMATION (LSA) 

Using H46|l in Ea. H20() allows the integrals over y° and K'^ to be evaluated as: 

S{x,K,q) = I 0^So{x;K°,K')e^^-^-'^^y^° J d^y e-'K'-''vI/(,;)(x + y/2)vl/(,;)*(x - y/2), (56) 
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FIG. 1: Plot of g(b), showing cutoff around b = 1/T, with T = l/m,r and b in fm. 



and using (|49|) allows the integrals over and K^' to be evaluated yielding a four-dimensional integral: 



S{x,K, q) 



;SoiT, ri)e 



iq t—iqiz 



J d%'B,{h, b')<;)(b + b72)^.(;)*(b - h'/2), 



(57) 



B^{h, b') = J d?K'j, B^{h, K'j,) exp [-iK^ • b'] . 



The result (|57() still requires the evaluation of an 6-dimcnsional integral (over r, ?7,b',Ky) to obtain 
the correlation function. We search for simplifications. The integral (|57|1 simplifies if we ignore the effects of 
transverse flow rapidity. So to gain insight, let's set 



Vf = 0, 

consider a fixed value of 77, and take one of the terms in the series H37|l . Then 



5(b'2) = 2 / d^Kj_ Mj_ exp 



—M± cosh 77 



T 



exp [— iK_L • b'] 



(58) 

(59) 
(60) 



Fig. fflshows g{b) for the case T/ cosh 77 ~ m^. It is clear from H6Q(I that the quantity T controls the range of 
allowed values of K±, so that the extent of b' is of order 1/T w 1 fm for T = 200 MeV. This is much smaller than 
the size of the presumed fireball that controls the extent of b in (|57() . Thus it is natural to think of neglecting 
the terms involving ±b'/2 of 1)57(1 . This however, is too extreme an approximation because it would not lead 
to the formalism of Sect. 1 in the plane wave limit. Instead, we use the approximation 



This is exact in the plane wave limit, but its wider validity relies on the replacement 

^(-)(b + b72)V;(-)*(b - b72) g(b'2) « ^p(,~\h)^p^^-> {h) ^(b'^) exp (^K^ • b'). 



(61) 



(62) 



that requires that the size of the source be much larger than n/T, where n is the expansion order in the Bose- 
Einstein distribution 1(37(1 . Our sources typically have a diameter of about 25 fm, and T ~ 1 fm^^, so that n 
must be no bigger than about 25. We achieve numerical convergence with n < 10, so that the source is truly 
large enough for our approximation. We denote Eq. 1(62(1 to be the "Large Source Approximation" (LSA) and 
use it to immediately integrate over b' and to obtain a simpler version of 1(57(1 : 



i{uj2 — LJi )r cosh — 



,.r.i„h,5^^(l,^K)7^(,;)(b)V'^,;^*(b), 



(63) 
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that is obtained for any value of 77/. The expansion ()37|l gives 
S^(b,K) = 2^exp(- 



n=l 



-)MTpib) 



E°° . -Mt cosh 7? cosh jyt (6) , sinh (&) cos </> 
exp( ) cxp( )MTp[b), 



n=l 

SO that Eq. H63() becomes 



Tn 



Tn 



(64) 



B„(b,K) EE exp(^^/T„)exp(/-Csinh77t(6) cos0/T„) Mt/9(6) 

OC 

5'(a;,Pj) =5o(T,7y,y)^cxp(-7„(&)coshry)B„(b,p,;) -^^y^b) 

ri=l 

Mrcoshijtib) 

7n(0) = . 



(66) 
(67) 

(68) 



The number of terms in the summation over 77 required to achieve an accurate result depends on the values of 
Ptt, K. If desired the quantities T, p^ can be treated as functions of b |l7j in Eq. (|66() . rather than as constants, 
as we have done here. 



VIII. CORRELATION FUNCTION 

The specific form of the wave function enables us to express the correlation function as C((/iZ + q, K): 

C(,.Z + q, K) . 1 + r,J{f'''['if:f , (69) 

Pi,2 = K ± (7i/2z ± q/2. (70) 

The dependence on occurs only in the numerator. 

The remaining task is to perform the integral over d'^x. We use two separate approaches. The first 
[Tol | involves expanding S{x,K), Eq. H66|l . as double power series in and r, keeping all terms up to second 
order. The second involves exact numerical integration. The two methods yield nearly identical results, with 
the second being more accurate and taking only slightly more computer time. We present both methods here. 
First, we take q to be small and make expansions Sec. IVIII Al This formalism is applied to obtain numerical 
results for radii and spectra in Sects. IIX AlIXEl and Sect. IIX Fl The formalism to compute the correlation 
functions without the use of expansion is contained in Sec. IVIII Bl and this formalism is applied to compute 
correlation functions in Sec. IIX Gl 

A. Evaluation of correlation function by expansion 

Now we make the above mentioned expansion keeping terms to order g£, (and anticipate the integration 
over 77, r). The term linear in is an odd function of rj so that it vanishes when the integration over 77 is carried 
out. Use a;2 — = —qoP, (i = Kt/Mt- Thus the expansion of Eq. is 

^(a;,^) =5o(T,77,y)(l-ig<,/3rcosh77- ig^^V^ cosh^ 77 - ig^T^ sinh^ 77) 



oc 

X 

n=l 



^exp(-7„(6)coshr;)i?„(b,KT)7A(,;nb)^^;'*(b), (71) 
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The effects of the Y dependence of So{T,rj,Y) that appears for ^ is the same for the numerators 
and denominators of the correlation functions and vanish. Hence we do not keep this exphcit dependence in 
the intermediate steps in the following calculations. 

The next step in evaluating H69() is to integrate over all r, 77 using the measure drfrdT. The first integral 
to appear arises from the factor of unity appearing inside the parenthesis of H71() . It is 



27rA7 



drj cosh rj I rdr cxp{' 



^ )cxp(- ^^ /""^ ) exp(-7n(6) coshT?) 



2(Ar;) 



2(Ar) 



« 2Toexp(^-2)/vi(7„ + ^-2) = To/o(Cn)exp(^--2), /o(Cn) = 2/^1 (61), 6i = In + 
The approximation involves the replacement 



(72) 



(73) 



Similar replacements are made below. The dominant contributions to the integral involve small values of 
ry, so the approximation is expected to be very good. The appendix shows that the error involved is less than 
about 2%. 

We proceed to use the same replacement to evaluate the remaining integrals and find 



1 



/27rA7 



cosh r]dri / rdr exp{- 



1 



I exp(- 



2(At)2 



) exp(— 7„(6) cosh 77)(—igo/3)T cosh ?7 



{-iqoP){T^ + (Ar)^) exp( — )/2(60, ./2(60 = i^o(Cn) + ^'2(6.) = 2(A'o(en) + 



V27rAT J- 



/OO 
rdr Gxp 
-00 

1 



(74) 



7]^ 



exp(- ° ) exp(-7„(6) cosh77)(qo/3)^T^ cosh^ r/ 



(3roAr^ +To^)(go/3)'cxp( — )/3(e„),/3(en) = 2 

^2 



2(A7y)2' 2(At)2 

^2(6.) 



2( A-o(en)/en + ^i(en)(i + 7j) 

Sn 



/i = / d77COsh?7 / TdTCxp{' 

V27rAT J_oo ^-00 



2(A^'' ^^^^ — 2(At)2 ^ '^^P("^"(^) cosh7])T2 sinh^ 77 



1 



ro(3Ar^ + t^) exp(^)/i(60, /i(60 = 2(^^2(6^)/^^) = 2/^0(60/^. + 4ifi(^„)/e^. 



(75) 



(76) 



Then 



j d'^xS{x, K, q) = To exp(^^) 



$12 - iqoPiro + At^/to)F2{Kt) ~ -qlP^SAr^ + T^)F:i{KT) - ql/2Fi{KT) 



where 



00 ^ 

^i2 = Y. rfVo(e„)i?„(b,KT)V^;)(b)4;)*(b) 

00 „ 

F2{Kt)^Y1 / d%h{in)Bn{h,\^T)\^^^l{h)\'' 

n=l 
00 „ 

F^{Kt) = Y. / d26/3(e„)Sn(b,KT)K'(fj(b)p 

n=l 
00 „ 

i^i(i^T) = E / rf'^'/i(5n)i3n(b,KT)|7Air)(b)p 

n=l 

The denominator is obtained by evaluating the emission function H()8|) and the following functions enter: 

Fo{Kt)=Y^ J dHfo{^n)Bnih,KTM^-^{h)\^ 



(77) 

(78) 
(79) 
(80) 
(81) 

(82) 



ri=l • 
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CO „ 

^" = E/ dHmn)B^{h,KTM^-'>ih)\\ (83) 

n=l 

With these definitions, one may show: 

C(q + 9^2, Kt) = 1 - - qlfS^Ar' + (84) 

'3?ll5>22 

Rf = (3Ar2 + T2)i^i(i^7,)/Fo(i^T) (85) 



At = (3Ar2 + t^)F3{Kt)/ FoiKT) - (tq + Ar^Tof 



F2[Kt) ' 



Fo{Kt) 



(86) 



The spectra are given by Eq. Q 0| 



so that in general 

dN 1 



(27r)) / / d''xS{x,p), 



in which the azimuthal symmetry of the angular distribution is used. 

The STAR detector receives pions for values of Y between ±0.5 and presents its results in terms of 
' '^^ )|r|<o.5 which is given by 



dN \ 



2TTM±dM±dY I 



dY I d''xS{x,K). (89) 

-0.5 



Numerical studies showed us that the average over Y is extremely well approximated by simply replacing Y by 
0. This requires 5ri > 1/2, and we use Srj « 1. Thus we find 

The evaluation proceeds by reducing the two dimensional integrals of Eqs. 1)78183(1 to those of one dimen- 
sion by an analytic evaluation of the angular integrals. We need the partial wave expansions 



^^;^*(b) = /o(pi,6) + 2 /m(pi,fc)(-0™ cosm0i (91) 

rn— l.oo 

i^i-\h) = r,{p2.h) + 2 /,„*(p2,6)(z)" cosm</)2, (92) 



where cos 4)i = pi • b. 

We encounter integrals of the form 



m— l,oc 



27r 

z cos ( 



Amn{z)= I dcf) 6^''°'"'' COS m4)i COS n(t)2- (93) 

JO 

For q II K, q = qo, 0i = 02 = 4>^ so that we define 

A^n{\\,z)= d0 e^^^-^cosm-^cosn^ (94) 



= TT {Irn+n{z) + Im-niz)) , (95) 

where Ini±n{z) are modified Bessel functions: 

Iniz) ^ {-irj„{iz), (96) 
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for real z, and In{z) = An integral representation is 



1 /■^'^ 



For q • K = 0, q = Qs, we define 



K 



cos a 



sm a 



g/2 



Then 



A,„„(±, z) = TT (/,„+„ (z) COS (m - n)a + /,„_„(z) cos (m + n)a) . 



(97) 



(98) 



(99) 



Note that for the denominators, one gets integrals in which the two angles of (|93|l are the same. Then the 
expression (|95|l is to be used. 

The use of (EU and (jHUl in (jHSjl yields 



$i2(i|,±) = M^{K) \ bdbpib)e 



-M±(K) coshijt{b)/T 



X 



e™e„/™(pi, 6)/:(p2, 6)(*)"-'"A™„((|l, ±), ^ sinhryt(6)) 



_m,n— 
Co = ^n>0 = 2. 



T 



The notation (||, ±) denotes either of the two possibilities q || K, q ± K. Similarly, 

/•OO 

$,,(||,±) = M^(p,) / 5d6p(6)e-^^-(f-)™^h''*W/^ 
Jo 

OO 

e™e„/™(p.,fe)/,:(p.,^)(*)"""^m4(ll,^),^sinhr;t(&)/r) 



.m.n— 



(100) 
(101) 



(102) 



B. Exact numerical evaluation of correlation function 

We present an alternate calculational technique that gives the correlation function for any value of q. To 
do this, start by recalling Eqs. (|69|) and H63() and note the appearance of the integral 



1 f°° 

= / (i77Cosh77exp(— 7„(6) coshyy) exp(— 

TtAt J-oo 

Trl-T rirnf ~ ''"o) \ iT((t^2 -t^l ) cosh Tj-igi. sinh?)) 

rdrexp^ ^^^^^^ )^ 

/OO 
dj] coshf] cxp{—^n{^) cosh 7/ 4- iaro) exp(- 
-oo 



r)(ro + iaAr^) exp(-a2AT72), 



a = {uj2 — uJi) cosh 77 — (7l sinh 77. 



The procedure of this section is to evaluate the integral over 77 numerically. 

Let's set up the full calculation. Integrating over r and using the result (|103|l yields 

n=l 

f°° IJ^ 

/ d77Cosh77exp(— 7„(6) cosh 7/ + iarg) exp( — - — r7T)(To + iaAr'^) exp(— At^ /2), 

J~oo 2(A77)2 



(103) 
(104) 



(105) 
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There are some remarks to be made here: there is no need to make the approximation of Eq. (|144|1 . and 
the present expression is actuaUy much more compact than H77|l . Possible cross terms j4r)j involving qolL are 
small for the experiments with y w that wc analyze, so we neglect the cross terms and take either (?l = 
or 90 = 0. Then the integrands have terms either even or odd in ry. The odd terms cancel. So define (with 



/(7„(fe), Atj,(7L, Ary, At,to) = 

/ (i?7Coshr/exp(— 7„(6) cosh?] + ickTo) exp(— — — — —)(to + ial^r )exp(— a Ar 12). 
J-oo 2{AT]y 



(106) 



We use various specific values of the arguments of the function / to compute the different observables. According 
to Eqs. (|25I26|I a radius Ri can be computed using qi ^ 0, qj^i = 0. Thus to compute Rq we take qL,s = 0, so 
that 



J(7, Au, 0, A?7, Ar, tq) = Io{l, Alu, At;, Ar, tq) 



2 / o?77cosh77exp(— 7C0sh77) exp(— 
Jo 

z = Alu cosh 77. 
To compute Rl we take Auj = 0, so that 



2(Ar7)2 



' (ro + izAr^) exp(-z2 ArV2) 



(107) 



7(7, 0, gi, At/, Ar, ro) = 7^(7, gi, Ary, Ar, ro) = 

2 / d77COsh?7exp(-7Cosh7?)exp(-— — — r)(rocos(j/ro) - yAr^ sin(?/ro)) exp(--?/^Ar^/2) 
Jo 2{Ariy 

y = -qL sinh?7 



(108) 



In computing the spectra, or the denominator of the correlation function we have Auj ~ 0,qL — 0. Then 
WC use: 



/(7, 0, 0, A?7, Ar, ro) = 2ro J d-q cosh 77 exp(-7 cosh 77) exp(- ^^^^^2 ) ■ 



Then the correlation function is given by 



C(q + gLZ,K) = 1 



IX12I 
X11X22 



X12 = j d'^xS{x,K,q) 

= Y. I rf'foi?n(b,KT)/(7n(6),Aa;,(ZL,A77,Ar,ro)7/'(,;)(b)7/-(,;)*(b) 

n=l 

00 „ 

X^^ = Y. j rf'^5n(b,p.)/(7n(6),0,0,A77,Ar,ro)|7/.(,-)(b) 



(109) 
(110) 

(111) 

(112) 



For very small values of q this correlation function reduces to the one obtained in second order, 1)84(1 . There 
are two changes. If q is not small then one obtains the correct correlation function. A second change is that 
the approximation (| 144(1 is not used. The differences between using the procedure of this subsection and that 
of the previous subsection are negligible (< 2% and indistinguishable) for the cases we have studied. This 
provides a further verification of the approximations used in Sect. I(VIII A(l . but the present formalism avoids 
those approximations. 



IX. APPLICATIONS 



The description of the formalism is essentially complete. The plane wave emission function 5o is defined 
in Sect. II and its symmetry discussed in Sect. IV. The distorted wave emission function function S is defined 
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TABLE I: Best fit parameters (F193) used in the calculations shown, with variances; fits use Eqns. (123123. Bold-face 
indicates parameters not varied in the fitting procedure (see text). 



r 


Vf 


At 


Rws 




Wo 


W2 


To 


A77 




(MeV) 




ifm/c) 


ifm) 


ifm) 


ifm-') 




ifm/c) 




(MeV) 


193.00 


1.507 


2.382 


11.773 


0.909 


0.139 


0.847 +i0.120 


8.60 


0.991 


0.000 139.57 




±0.025 


±0.07 


±0.06 


±0.015 


±0.046 


±0.014 ±0.002 


±0.10 


±0.032 





in Sect. Ill, and its evaluation elaborated in Sect. VII. The optical potential is defined in Sec. V, and its use in 
a wave equation to obtain the distorted wave is discussed in Sec. VI. The correlation function is evaluated in 
Sect. VIII. 

Our technique may be compared with the Buda-Lund model, an efficient representation of the data|25). 
in which the temperature and fugacity are taken as position-dependent functions appearing in a Boltzmann 
distribution. The effects of our optical potential could provide an explanation of some of those deduced depen- 
dencies. 

We are now ready to confront the data. 

A. DWEF Fits to Central Au+Au Collisions 

Table 1 gives the parameters of our best fit to the STAR data for Au+Au central collisions at 200 GeV, 
which we will refer to as "F193". For this fit to the data, as shown in Figs. 13 and ^ the is 56.45, and 
the per degree of freedom is 2.45. Table I also gives the estimated variances of those fit parameters that 
were varied, as calculated by determining the parameter variation required to increase the value by one unit. 
Correlations between different parameters are not considered, but, as will be discussed below, more than one 
set of parameters can produce a quality fit to the data. 

To illustrate the importance of the various effects of the DWEF model for computing the radii Ro,Rs 
and the spectrum, we have done calculations using the F193 fit parameters with various effects switched on 
separately. This is shown in Fig. [51 The curves labeled DWEF show the full calculation. Those labeled PWEF 
are computed using plane waves, i.e., the optical potential (^0,2) and flow (rjF) are set to zero. The curves 
labeled Re{wo) and Re{w2) use only the real constant or momentum-dependent parts of the optical potential, 
respectively, set to the F193 values, with no flow. The curves labeled Im{w2) use only the imaginary momentum- 
dependent part of the optical potential set to the F193 value, with no flow. Finally, the curves labeled flow 
uses the F193 value of -qp, with no optical potential. All other parameters in these calculations are set to the 
F193 values of Table 1. This study indicates that both flow and the optical potential modify the HBT radii, 
but only the momentum dependent parts of the optical potential (1^2) affect the spectrum. 

Figs. 121 and ^ show DWEF calculations using F193, as compared with the STAR data |2^I3 used in the 
fit. We see that both the HBT radii and the spectrum (including the spectrum normalization) are reproduce 
very well indeed. However, some clarification is needed here on the issue of resonance pions. Our model 
describes only those "direct" pions that are emitted directly within the fireball and that participate fully in the 
HBT correlation. It does not predict the "halo" fraction of pions originating from resonance decays later in the 
process that will not participate in the HBT correlation (or that would make an unmeasurable "spike" in the 
correlation function near q ^ Q). These latter pions are present in the measured spectrum and must be removed 
before fitting the spectrum with a DWEF calculation. 

We have done this by using the HBT A parameter previously extracted as part of the HBT analysis, 
which we take to represent the probability that both correlated pions are direct rather than halo pions. We fit 
the extracted A parameter with a straight line in transverse momentum to obtain the function A(pt)- Then we 
remove the 12% weak decay correction of the published pion spectrum |34j and multiply it by \/'\{j>t), thereby 
producing an estimate of the momentum spectrum of direct pions only. This corrected spectrum, as indicated 
by the black small-dashed line in Fig. 01 is then used in the data fitting. Following the DWEF calculation, the 
predicted theoretical spectrum is "uncorrected" by dividing it by y/Mpr) and then reducing it by 12%, so that 
it can be compared directly with the published data. These predictions, now including halo pions, are the curves 
shown in the spectrum plots. We note that the effects of the two corrections to the spectrum tend cancel each 
other, and that their net effect on the fits is primarily an increase in temperature and fiow parameters. We note 
also that the Phobos Collaboration has published spectrum data in the lowp^ region|3^, but because the value 
of A in this region are not known, we cannot make similar corrections and therefore have not included these 
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FIG. 2: (Color online) Calculations of Ro, Rs, and spectrum, isolating potential and flow effects. DWEF (black 
solid) full calculation; PWEF (red dashed) plane wave calculation with no optical potential or flow; Re{w2) only (green 
dot-dot-dash) calculation with no flow and only real part of momentum dependent optical potential (1U2); Re{wo) 
only (blue dot-dash) calculation with no flow and only real constant optical potential (wo); Im{w2) only (solid violet) 
calculation with no flow and only imaginary part of momentum-dependent optical potential {W2); Flow only (orange 
dotted) calculation with flow but no optical potential. Non-zero parameters in all calculations have the F193 values of 
Table I. 



data in our current analysis. The DWEF fit parameters, particularly the temperature, flow, chemical potential, 
and Woods-Saxon diffuseness, are sensitive to the details of our treatment of resonances, thus introducing a 
procedure-dependent systematic error in DWEF parameter extraction. 

Our DWEF calculations predict the correlation function of two pions with a particular vector momentum 
difference q. There arc at least two ways of extracting HBT radii from such correlation functions. One is to 
use the quadratic form H24|) of the correlation function at very small magnitudes of q. The other is to use the 
Gaussian form H23|l at larger values of q that correspond to the falloff region of the correlation function. In our 
previous publication^^! we used the former method with q ~ Kt/'^Q, a technique which essentially extracts 
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FIG. 3: (Color online) HBT Radii Rs,Ro„ Ri and the ratio Ro/Rs; Data V (green) tv+tv+; A (red) tv-tv'. 

Curves: solid (black) full calculation; dotted (red) t?/ = (no flow); dashed (orange) Re[f/]=0 (no refraction); 
dashed (green) U—0 (no potential); dot-dashed (blue)— >, = (almost same as solid) double-dot-dashed (violet) 
substituting Boltzmann for Bose-Einstein thermal distribution. 




100 200 300 400 5Q0 600 700 

Kr (MeV/e) 

FIG. 4: (Color online) Pion momentum spectrum. Data Is^l : V (green) =J> tt"*"; A (red) tt^. Black dotted line is the 
'non- resonant pion" spectrum used in the fit. The curves are denoted as in Fig. El 
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FIG. 5: (Color online) Wave functions. The figures show the absolute square of the calculated wave functions times the 
density p{b). For each (row) value of Kt (100,250,600) MeV/c, the full calculation, the calculation including only the 
imaginary part of the optical potential and the eikonal approximation with the same absorption are shown horizontally. 
The out direction is parallel to the lower right edge of each picture. 

HBT radii from the calculated curvature of the correlation function near q = 0. In the present work we use 
the Gaussian form with q — 0.15 fm^^, a momentum difference that evaluates the correlation function near 
its half-maximum point, a procedure that bears more resemblance to the extraction of radii from experimental 
data with Gaussian fitting. We have found that in the medium and high momentum regions the two procedures 
gives very similar descriptions of the data and extract similar radii. However, we have found that for values of 
Kt < 175 MeV/c, the two methods diverge, and that the Gaussian form gives more reliable radii. 

B. Wave Function Plots 

We can understand more about how HBT measurements in various momentum regions probe the system 
under study by examining the computed wave functions from the DWEF calculation. Such wave functions 
calculated with the F193 fit are shown in Fig. |31 The complete DWEF wave functions (left column) differ 
significantly from the semi-classical eikonal approximation calculated with the same absorption (right column) 
at each of the energies we consider. The importance of the real part of the optical potential is illustrated 
by comparison with the center column, showing wave functions calculated with the real part of the optical 
potential set to zero. The maxima and minima of these plots result from interference effects caused by the 
optical potential as incorporated by solving the quantum mechanical wave equation. Their spacing gives a 
qualitative indication of the pion wavelength in the medium. The differences between the quantum calculations 
and the semi-classical eikonal model grow smaller as the value of Kt increases. We note that while the high 
momentum wave functions sample only a "bright ring" at the edge of the fireball, the low momentum wave 
functions is also non-vanishing in most parts of the fireball. 

C. Temperature Ambiguities 

In our previous publication |lOl| the temperature parameter used in the pion emission function was T = 222 
MeV. This is an uncomfortably large value for the temperature of the medium, high enough that pions would 
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TABLE II: Nine Fit-Parameters Sets: These fits were obtained by fixing a range of values for the temperature T and 
fitting the central Au+Au 200 GeV STAR data. Bold-face indicates a parameter fixed for the search. 
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be expected to "melt" in such an environment and lose their identity. Therefore, we explored (not shown) 
the extent to which the DWEF model requires such large values of temperature and showed that there is 
a parameter ambiguity involving the parameters T and awS: such that a fit of reasonable quality could be 
obtained for temperature values over a fairly broad range. 

Therefore, it is desirable to choose a temperature appropriate to the medium. If one takes quite literally 
the expectation that the DWEF model describes the initial emission of pions and that the first pions are 
produced directly in the strongly interacting quark-gluon plasma as it make a phase transition to a hadronic 
phase, then the emission function should have the temperature of the QGP transition. Similarly, if the pions 
are produced as massless objects due to chiral symmetry restoration in the medium, then the chemical potential 
should equal the free pion mass. Recent lattice gauge calculations reported at Quark Matter 2005^ give the 
critical QGP transition temperature to be 193 MeV. Therefore, we adopt (Table I) T=193 MeV and ^7r=139.57 
MeV and search for a new fit to the STAR ^/5jvjv=200 GeV Au+Au data. To our surprise, the fit we obtain 
(see Figs. 3 and 4) with parameters fixed at these values is the best we have found, with an overall of 56.45 
and a Per degree of freedom of 2.45. Further, subsequent searches in which the fixed temperature was set 
to other values between 173 MeV and 220 MeV show that there is a definite minimum in at just the value 
of temperature given by the lattice gauge calculation. Table II shows nine different set of fit parameters, all 
of which give good fits to the central STAR ^/SNN~200 GeV Au+Au data. In Table II, the fixed parameters 
are indicated in bold face. We see that F193, the parameter set of Table I, gives the best fit, and that for 
temperatures greater than 193 MeV it is necessary to use a non-zero value of e, invoking a Kisslinger-type wave 
equation Eq. H45|) . to obtain a quality fit. However, fits F193a and F193b indicate that searching on e as part 
of a search at T = 193 MeV gives values near zero and did not improve the fit. 

Fig. |B1 shows per degree of freedom as a function of temperature for the nine fits listed in Table II. 
The solid line is a parabolic fit to the points. We do not assert that there is physics driving the preference of 
the STAR data for the temperature predicted by lattice gauge calculations, but we believe this is more than a 
coincidence. 

Fig. Ul shows the ratio of the parameters of Table II to the F193 parameters, and indicates the range of 
variations and the correlations of parameters. Note in particular the correlation between T, aws, wo, Ar, and 

VF- 

Fig. |S1 shows the predictions of all nine fits as compared with the STAR data. As can be seen, there are 
no striking differences. In order to make the differences more visible. Fig. |51 shows the ratio of the predictions 
of the eight other fits to those of F193. as c an be seen the differences in the radius predictions are less than 
1%, while those of the spectrum predictions are as large as 10%. 

However, we note that in the momentum region below 50 MeV/c, where no HBT or spectrum data is 
available, the different fits make different predictions. Fig. 1101 shows the predictions of the fits of Table II in 
this region. We note that the "wiggles" arise mainly for the extreme fits, and in particular that these variations 
are relatively small for the F193 fit that provided the best fit to the higher momentum data. 
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FIG. 6: (Color online) Dependence of chi-squared per degree of freedom on temperature. Symbols correspond sequentially 
to the fits listed in Table II: circle (red), plus (orange), triangle (green), diamond (chartreuse), 5-star (aqua), 6-star (light 
blue), 7-star (dark blue), square (violet), and del (magenta). Black line is a parabolic fit to the points, provided to guide 
the eye. 



D. Meaning of optical potential parameters 

Table I shows that wq = 0.14 and that W2 is 0.85 +0.12 i. Let us try to understand these values using 
the ideas of Sect. V. First consider the real parts. Comparing Eq. H45|l and Eq. H44() shows that 

TJV^ 1 — li"^ 

wo^^-ml{T), Re(u;2) = (113) 

which allows us to extract the values: 

ul = 0.54, = 0.74, u^m^iT) = 0.65 fm~\ (114) 

These values are comparable to the estimates of 13(1. 

We also compare with the results of Shurvak|37j for pions propagating in hot dilute matter. His results 
for T = 200 MeV (the stated limit of validity of the calculation) are most relevant for our case. His Fig. 1 
shows results for V ~ J7/(277i^) as a function of pion momentum. The imaginary potential is approximately 
proportional to for p < AOOMeV/c. For p = 1 fm^^ Shuryak obtains Im{U) « 0.1 fm^^, which is close to our 
value of 0.12 fm^^. Our real potential is about seven times larger in magnitude than our imaginary potential, 
while Shuryak finds that the strengths are comparable. Thus our real potential is much stronger than expected 
from a dilute gas approximation and consistent (via the formalism of |30j ) with the occurance of a chiral phase 
transition. 
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FIG. 7: (Color online) Ratio of other fit parameters to those of F193 for the nine fits listed in Table II. Symbols 
correspond sequentially to the fits listed in Table II: circle (red), plus (orange), triangle (green), diamond (chartreuse), 
5-star (aqua), 6-star (light blue), 7-star (dark blue), square (violet), and del (magenta). 




100 200 300 400 500 60O 700 
K,(MeV/c> 



100 200 300 m 500 600 790 
K, (MeV/c) 



FIG. 8: (Color online) Superposition of nine fits with STAR data, showing that there is relatively little difference 
between the predictions. Data points: V (green) => tt"''; A (red) tt~ . 



E. Non-central Au+Au 



Our analysis so far has focused on the central (0-5%) STAR v^5Arjv=200 GcV Au+Au data. However, 
the STAR collaboration performed measurements of pfon correlations and spectra at VSww=200 GeV Au+Au 
function of centrality. 

For non-central events, our optical potential would depend on the direction of b as well as its magnitude. 
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FIG. 9: (Color online) Ratio of other predictions to F193. Colors correspond sequentially to the fits listed in Table II: 
red, orange, green, chartreuse, aqua, light blue, dark blue, violet, and magenta. 
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FIG. 10: (Color online) Low-momentum predictions of the nine fits listed in Table II. Colors correspond sequentially to 
the fits listed in Table II: red, orange, green, chartreuse, aqua, light blue, dark blue, violet, and magenta. Note that the 
low temperature fits (red to aqua) have relatively little structure. 



The simple dependence on b was exploited heavily in previous sections to simplify the calculations, so, in 
principle, non-cental collisions are not a part of our model. However, we can make simplifying assumptions 
that can allow us to predict the observables for non-central collisions. In particular, we assume that a non- 
central collision resembles a central collision with the same number of participants. This assumption allows 
us to extrapolate our results to systems that do not have perfect centrality by using participant scaling. In 
particular, we take the space-time parameters RwStO,ws^ a-nd tq to scale as the centrality-dependent number 

1/3 

of participant particles to the one-third power: N^i^^^. The values of A^part arc taken from Glauber- model 
calculations |3^. For the Au-|-Au system, the value of At is kept at the value of F193, because this is a dynamic 
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TABLE III: Au+Au scaled space-time parameters from fit Ff93 
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TABLE IV: Cu+Cu scaled space-time parameters from fit Ff93 
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To 


Ar 
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ifm/c) 


ifm/c) 
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98.38 
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1.629 
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4.613 
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quantity describing the proper-time duration during which pions are emitted in the collision. Table III gives the 
parameters R\ys,'''Ot'^ws 3'iid Ar vs. centrality, as scaled from the F193 fit of Table I. Parameters not listed 
in the table are the same as those of Table I. 

In Fig. Illl wc sec that the participant-scaled predictions work fairly well for the first few centrality bins, 
but agreement diminishes for the more peripheral collisions. We do not expect this simple scaling procedure to 
be accurate very far from the purely central collision assumption of the model, so this loss of predictive power 
is to be expected. 

In particular, looking for differences between peripheral and central collisions is a standard way to test 
for the possible existence of QGP effects. This means that the strength of the optical potential parameters of 
our model should change as the collision becomes peripheral. We have kept the potential depth parameters 
constant with centrality because we have no way to predict their centrality dependence. Fig. 1111 shows an 
excellent reproduction of R^ with centrality, a fairly good reproduction of Rq, but only a fair description of 
Rs- We take this to be consistent with the disappearance of a chiral phase transition in the more peripheral 
collisions. 



F. Predictions for Cu-|-Cu collisions 

The STAR collaboration has also measured (but not yet published) radii and spectra vs. centrality for 
Cu-I-Cu collisions at ^/Snn=2Q0 GeV, so it is of interest to make predictions for this system. Here again use 
participant scaling, but have also assumed that the emission duration parameter Ar scales as A^^^, where A 
is that atomic mass number of the colliding nuclei. Table IV gives the scaled space-time parameters used to 
predict the HBT radii for the Cu-|-Cu system. Parameters not present in the table are the same as those of 
Table I. 



G. Correlation functions and the Gaussian approximation 

We apply the present formalism to obtain correlation function C (with A = 1), using the parameters of 
Table I, for Kt = 158, 316 MeV/c. The results are shown in Figs. I13|T^ We see that the correlation functions 
C{Kt, q) are fairly well represented by Gaussians, for the relevant range of small values of with widths that 
are approximately independent of Kt- For a typical radius of about 7 fm, and q = 0.15 fm^^, (where data are 
measured) q^R^ w 1. Therefore, using the approximation (|24|) is not very accurate. On the other hand, the 
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FIG. 11: (Color online) Centrality dependence of HBT radii in Au+Au collisions. Curves are predictions from 
participant-scaled space-time parameters; data points are Au+Au at y^sjvjv = 200 GeV from . Symbols and 
colors: 0-5% circle (red), 5-10% plus (yellow), 10-20% triangle (green), 20-30% diamond (aqua), 30-50% 5-star (blue), 
50-80% 6-star (indigo). 



correlation functions are close to Gaussian in shape in this region, suggesting that approximation H23|l is more 
appropriate. A detailed look at the ratio of computed correlation functions (C— 1) to its Gaussian fit in Fig. 1141 
shows that the Gaussian curves represent the correlation functions fairly well in the region < q < 0.22 fm"-'^ 
where C — 1 is large and radii are extracted, but that the correlation functions are systematically larger than 
the Gaussian fit for Ro,l and smaller for Rs in the "tail" region {q > 0.22 fm). 



H. Possible Extensions of the DWEF Model 



The DWEF model presented here uses the empirical "hydrodynamics-inspired" emission function of 
Eq. (|?7|) . However, we note that the apphcation of distorted waves to an emission function is more generally 
applicable, and that the formalism we present here can be applied to any emission function that has the same 
symmetry properties as Eq. I|27|l . One extension of the model would be to calculate the emission function as a 
multi-dimensional numerical table directly from hydrodynamics and use this with the DWEF model to calculate 
spectra and radii. Wc also note that, while it has not been not implemented here, the DWEF model can allow 
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FIG. 12: (Color online) HBT Radii; Central Cu+Cu = 200 GeV) predictions: Au+Au F193 fit - black solid; 
Cu+Cu Centrality 0-10% - red solid, 10-20% - yellow dot, 20-30% - green short-dash, 30-40% - aqua long-dash, 40-50% 
- blue dot-dash, 50-60% - violet dot-dot-dash. 
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FIG. 13; (Color online) Central Au+Au (^=200 GeV). Correlation functions (A't=158 , 316 MeV/c) for the out, side 
and longitudinal cases. The solid (black) curves are the full correlation functions. The broken curves are Gaussian fits: 
dash for out, dot for side and dot-dash for longitudinal. 




FIG. 14: (Color Online) Central Au+Au (^5=200 GeV). Ratios of correlation functions {Kt=158, 316 MeV/c) to 
Gaussian fits for the out, side and longitudinal cases. The dashed (red) curves are obtained with q in the out direction, 
the dotted (green) curves with q in the side direction, and the dot-dashed (blue) curves with q in the longitudinal 
direction. 
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FIG. 15: Eikonal calculation of pion wave function. A pion produced at a position b propagates along a straight line 
path 1 until it hits the edge of the medium, denoted by R. The angle between the dashed line and the horizontal line 1 
is denoted in the text as 9. The out case with = is shown. 



the temperature and chemical potential to depend on 6, as is done in the Buda-Lund model |25j. Prehminary 
investigations of this extension of the model, however, did not lead to an improved description of the data. 



X. THE EIKONAL APPROXIMATION 

Several previous calculations [iJl-E^I effects of opacity have used the eikonal approximation to 

solve the wave equation (|50|l . Here we explain the nature of this approximation and discuss its weaknesses and 
(fewer) strengths when applied to the current situation. 

The basic idea is that if the momentum is large one may say approximately that the wave propagate in 
a given direction (here the out direction, which is taken as along the x axis). Then one assumes a solution of 
the form ■!/)'^~)*(b) = e*P^<i>(b) that is inserted into the wave equation. Taking the Laplacian of the approximate 
wave function gives a term proportional to that is canceled, a term proportional to p that is kept, and another 
term that is ignored. Then one finds 



exp 



2p 



U{x',y)dx' 



(115) 



for propagation in the x (out) direction. The corrections to this solution are of order -^^jj^ — times 
the terms that are kept. The first correction can be large in the surface region in which U varies greatly and 
the second term can be large in the interior region in which U reaches its full value. We are concerned with 
pions of momentum ranging from 30 to 600 MeV/c, so that the eikonal approximation be expected to be poor. 
However, the ease of application, and its wide use makes it worthwhile for us to assess the use of Eq. ()115|l . We 
consider a purely imaginary potential first and then use a general complex potential. 



A. Strong Absorption at High K — Purely Imaginary potential 



In the impulse approximation U = —A'njp where / is the projectile-target scattering amplitude and p is 
the density of scatterers. The optical theorem relates the imaginary part of / to the total cross section a so 
that Im[U] = ^pap. If we keep only the imaginary potential ~i/{2p)U = \/2ap and for a constant density the 
intensity of the wave falls as e~^'^'"'p with the mean free path, Amfp = More generally the wave function for 
a purely imaginary optical potential is given by 

'^l>^-\v^M = e-*P'-''e~''/2^»fp (116) 

where i = 1,2 for the two wave functions and li{h, K) is the direct line path length (parallel to the direction of 
Pi from the emission point b to the edge of the medium. For a purely imaginary optical potential U = —iKap 
and Anifp is the resulting mean free path. 
We have, sec Fig. 1151 



R = li + b, 

= ll + + 2bh cos {e - Os). 



(117) 
(118) 
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FIG. 16: Radii Ro and Rs for increasing values of Im['W2]. All other optical potential parameters and the flow parameters 
are set to 0. The numbers in the flgure refer to the range of values. The dashed curve shows the result of Heiselberg & 
Vischer. 



9 is angle between b and x-axis (direction of K) 9s is angle between pi_2 and K. This was called a in previous 
sections. The angle between b and li is 9 — 9s. Solving wc find 



/ = -fecos(6' - 9s) + y/{hcos{9 -9s)Y +R'^ -h^. 
For the R-out case 9s — 0, and 

U = l = - 



(119) 



(120) 



In this approximation the correlation function minus one is given by the absolute square of the ratio of integrals: 
particular, 



^||.i(<?=0)- 



/||(g) = y d^fee^i-^e-'/^-fp, q j] X 



A 



1 + iq\ 



(121) 
(122) 



H Amfp << R we neglect the second term. Since we are interested in radii, we expand the remaining term in 
powers of the exponential (keeping Amfp << R): 



A 



mfp 



igAnifp 



{2R - q^R^/2{2R - 2R/3) + iqTtR/2) 



G\\{q) 



g'i?'(2/3-7r7l6) 



i?^ = i?2(2/3 - 7rVl6) = 0.0498i?2; = R/AAS 



(123) 

(124) 
(125) 



In the plane wave approximation (Eq. \122\ with A,nfp oo) one would find R^ = Rj VS- The striking result 
of Heiselberg & Vischer [ijl indicated that the measured radius should be 40% smaller than the radius obtained 
in plane wave approximation. This result is confirmed for the highest value of Kt = 3 fm""'^ k, 600 MeV/c by 
the calculations shown above in Fig. |21 The value of i?o is reduced by approximately 40% by the influence of the 
imaginary optical potential. Fig. 1161 shows the effects of increasing the imaginary potential (by varying /m[w2] 
from 0.1 to 0.5 in steps of 0.1). The computed value of Rq does not change when lm\w'A is large enough and 
Kt is high enough. Small deviations between our results for highest Kt and the result |4y can be attributed 
to the non-zero value of the diffuseness aws ■ 

In the model of the present section, the correlation function and pion intensity would each be proportional 
to Amfp, and a very small value would yield a very small pionic spectrum. More generally, the extraction of the 
chemical potential from the pionic spectrum depends on the mean free path parameter Amfp- 
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Now let us do the R-side case. Then 



cos6'^ = ]sm0s = ±— (126) 



2 

cos (0-^.)=cos^?(l-^)T^ sine. (127) 



We find that 



k-^~x + ^B? ~y^± 0{j^) + 0{j^) ^-x+ - y2. (128) 

The terms ±0{-^) cancel in computing the term h+h that enters in computing the present correlation function. 
Thus the correction is of order 1/K^. The eikonal approximation works only if KR ^ 1, so the corrections 
must be presumed to be small. Then 



/■R. r\/R''-y'- i_ /^i— 7 

I±{q)= dye'iy dxe^/^-'pe ^»£p^^ ^ (129) 

« A / dy{\ - g2j;2/2) ^ A,„fp(2i? - q^R^/i) (130) 

J-R 

C^{q)^2-q^Ry3 (131) 
Rl = i?V3, (132) 

Without distortion would be (q) ~ 2 — q^R^ /A, so in this case the strong absorption increases the radius by 
a factor of ■\/4/3. This result is qualitatively obeyed in our realistic solutions of the wave equation, see Fig. 1161 



B. Complex potential with non- vanishing real part 

If there is an attractive real potential the wave function of Ea. (|116|l becomes: 

ip^-^pi^h) = e-*P^-'"e-'i(^+'")/2^»fp (133) 

with a dimensionless, real and positive. If a = one obtains the purely absorptive model of the previous 
sub-section. Conversely, the limit of a purely real potential occurs when l/Amfp 0, a/Amfp I/Aq. We also 
have 

In the product V'piV'pa enters in computing the correlation function, so unlike the previous case of pure absorp- 
tion, a term of the form li — 12 enters. This difference is of order q/K compared to other terms, but its influence 
in computing radii must lead eventually to a term of order (■§^)^. The validity of the eikonal approximation 
depends on the ability to disregard such terms. Thus a valid eikonal approximation means that h = I2 so the 
factors of a cancel out in the product ■0*-~^(pi, b)-0*^~^(p2, b). Thus, if one assumes the eikonal approximation 
is valid at all values of Kt, one would find erroneously that the real potential never has an influence on the 
calculation of HBT radii. 

Unlike the effects of the imaginary potential, which are qualitatively captured by the eikonal approxima- 
tion (even if applied wrongly at low values of Kt), the effects of the real potential are completely lost. Thus 
the eikonal approximation can not be used for values of Kt such that the real potential contributes. As shown 
in Fig. O the real potential is important for all values of Kt less than 600 MeV/c. Conversely, for much larger 
values of Kt, for which the eikonal approximation does accurately reproduce the solution of the wave equation, 
the real potential will not play a role in determining radii. 



XI. OSCILLATIONS - A SIMPLE SQUARE WELL EXAMPLE 



Our numerical results indicate that the radii may have significant oscillations for small values of Kt = K. 
The purpose of this sub-section is to provide a simple example that also yields oscillating radii. Consider a 
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cylindrical source of radius R of infinite extent in the longitudinal direction. Suppose this leads to a real, 
attractive square well potential of radius R that is proportional to the square of the momentum, as motivated 
by Eq. (gSJ. Then ^SU^ becomes 

(-V^ - C/op2)^ = (6 < R) (135) 

-VV=pV^ {b>R), (136) 

with Uq > 0. This equation is easily solved using the partial wave expansion H53(l . To provide a simple analytic 
example we further specify to the case of the lowest partial wave, m = 0. In this case we find 



ij{h) = MpVI + Uo 6), (b < R) (137) 

This shows immediately that the effect of the interaction is to scale each momenta pi,P2 that appear in the 
wave functions by a factor ^/l + Uq. We define 



p = p^l + Uo. (138) 

Eq. H137|l is a valid approximation to the full wave function (for b < R) only if p^/1 + UqR <C 1. However, it is 
interesting to also consider larger values of p. 

To compute radii, recall the correlation function Ht)9|) . For simplicity we take rjf — 0, and consider fixed 
values of rj and t, with At = 0. This corresponds to evaluating 5*0 at its peak and neglecting the influence 
of time duration. In this case many factors in the numerator and denominator of Eq. (|70|l cancel. Then the 
correlation function is given by a simple expression that provides some insight: 

(139) 

(PbXPi)<Pr{P2) 

, I \ - ^ f^'uMT f u\ T r u\ PiJa{P]R)Ji{PiR)-PjJo{piR)Ji{P]R) 
4'r{Pi.P]) = I bdbJa{p,b)Ja{pjb) = -^y^ — — (140) 

(141) 



R'J, R{PI-P^) 



■iP^) ^ ^2 £ bdhjlip^b) = i {4{P:R) + Jl[p^R)) 



If we are concerned with the side radius, inside the well Pi = P2 = \J K"^ + 9^/4 for the arguments of Jo that 
enter. In that case, (j)ii(pi,p2) 4'r{pi,pi) = 4)r{pi) and the correlation function takes on the value of 2 and 
the side-radius vanishes. This is a specific consequence of the approximation of taking only m = 0-all directions 
of K are equivalent, so there is no infiuence of vectors q that are perpendicular to K. 

The out case is more interesting because the energies and magnitudes of momenta pi.2 = {K ± q/2) of 
the two pions are different. A non-zero radius is obtained by evaluating C{K,q) for very small values of q and 
using Eq. The quantity Ro{K)/R is displayed in Fig.lTTI 

We see that an oscillatory pattern similar to our realistic results, such as Fig. 13 for small values of Ktj 
emerges from the present calculation. This occurs simply because Bessel functions oscilate. Indeed the out- 
radius obtained without the influence of final state distortions (t/o = 0) also has oscillations, but with a different 
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frequency. Therefore another quantity of interest is the ratio of radii computed with Rq (K) and without Rq (K) 
the influence of distortions: 

displayed in Fig. El The oscillations seen in the right part of Fig. El demonstrate the significance of final 
state interactions that cause enhancements by factors of greater than six and suppressions by factors of 2! 
Furthermore, comparing the left and right hand sides of the figure shows the quantity Rq{K)/R also oscillates. 
This is caused by the sharp edge of the square well and because we limit ourselves to to = 0. 



XII. SUMMARY 



A complete formal treatment of the distorted waves treatment of HBT correlations is presented here. 
The need for incorporating the influence of an optical potential U and the resulting distorted wave emission 
function is explained. The partial wave formalism necessary to compute the pionic distorted waves and the 
resulting emission function is detailed. Two different methods with equivalent results to evaluate the necessary 
eight-dimensional integral are described. Chiral symmetry restricts the form of U [sol l3ll| at low energy and 
the necessary constraints are implemented here and in [iQ. An excellent description of the STAR Au+Au 
HBT and spectrum data is achieved for central collisions and the use of an average area formulation leads to 
a very good description of these observables for non-central collisions. We also use four different versions of 
geometrical scaling to predict the results of central Cu-|-Cu collisions. The ability to calculate the absolute 
magnitude of the spectrum as well as the radii Eq. H25|) .Ea. H26|l (which involve ratios of functions of emmission 
functions) required for the computation of radii is a principal advantage of our method. The Blast Wave Model 
is discussed in Sec. IIVI 

The above results are achieved using a temperature Tc fixed at the recently determined critical value 
of 193 MeV Such a value could present difficulties for conventional calculations of the spectra because 

chemical equilibrium analyses yield lower temperatures Tch = 174 MeV ^ large difference between Tc 

and Tch implies that the hadrons interact after the deconfinement transition occurs. This notion is entirely 
consistent with our treatment of pionic distortions which has as its fundamental assumption that pions interact 
in a hot dense medium before escaping to freedom. 

We find that the necessary real optical potential is so strongly attractive that the pion can be said to 
lose its mass inside the medium. That chiral symmetry seems to be restored is the conclusion of our earlier 
work[l3. The RHIC-HBT puzzle is therefore replaced by the need to investigate this restoration. 

Explicit evaluation of wave functions obtained by exact numerical solutions of the wave equation show 
some interesting features of the strong interaction and also display differences with the solutions obtained using 
the eikonal approximation. A critical discussion of the cikonal approximation as applied to computing HBT 
radii shows that its use causes the crucial influence of the real part of the optical potential to be entirely lost. 
The huge importance of the real part of the optical potential is explicitly illustrated through a simple 

There are many immediate applications of this formalism. In particular, a treatment of HBT data 
obtained at all energies is in progress and will be presented elsewhere. 
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Appendix 

Consider the integral appearing in Sect. IVIII Al 

/>oo 

/(/3,A?7)=/ d?7 exp (-?7V(2A772)) cosh (t;) exp (-/3 cosh 77), (143) 
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Ratio 




/3 



FIG. 18: The ratio of the approximate expression Eq. 114411 to the exact one of Eq. (11431 . The values of 7 range from 
.05 to 0.2, with ratio approaching unity as 7 approaches 0. 

that can be approximated using Eq. (|7!^ as 

/(/3,A77) = 2exp(^)i^i(/3+^). (144) 
Figure ITSl shows that the approximation is excellent. 
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